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ABSTRACT 

The Q/U Imaging ExperimcnT (QUIET) has observed the co smic microwave background (CMB ) 
at 43 and 95 GHz. The 43-GHz results have been published in IQUIET Collaboration et all (|2011[ ). 
and here we report the measurement of CMB polarization power spectra using the 95-GHz data. This 
data set comprises 5337 hours of observations recorded by an array of 84 polarized coherent receivers 
with a total array sensitivity of 87/iK-^/s. Four low-foreground fields were observed, covering a total 
of ~ 1000 square degrees with an effective angular resolution of 12. '8, allowing for constraints on 
primordial gravitational waves and high-signal-to-noise measurements of the E'-modes across three 
acoustic peaks. The data reduction was performed using two independent analysis pipelines, one 
based on a pseudo-Cf (PCL) cross-correlation approach, and the other on a maximum-likelihood 
(ML) approach. All data selection criteria and filters were modified until a predefined set of null 
tests had been satisfied before inspecting any non-null power spectrum. The results derived by the 
two pipelines are in good agreement. We characterize the EE, EB and BB power spectra between 
I = 25 and 975 and find that the EE spectrum is consistent with ACDM, while the BB power 
spectrum is consistent with zero. Based on these measurements, we constrain the tensor-to-scalar 
ratio to r = 1.1±2J (r < 2.8 at 95% C.L.) as derived by the ML pipeline, and r = l.2to% (r < 2.7 
at 95% C.L.) as derived by the PCL pipeline. In one of the fields, we find a correlation with the 
dust component of the Planck Sky Model, though the corresponding excess power is small compared 
to statistical errors. Finally, we derive limits on all known systematic errors, and demonstrate that 
these correspond to a tensor-to-scalar ratio smaller than r = 0.01, the lowest level yet reported in the 
literature. 

Subject headings: cosmic background radiation — Cosmology: observations — Gravitational waves — 
inflation — Polarization 
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1. INTRODUCTION 

The theory of inflation expl ains several well-obse rved 
properties of the Universe (e.g. iLiddle fe Lvthll2000l and 
references therein) : the lack of spatial curvature, the ab- 
sence of relic monopoles from a grand unified theory's 
broken symmetry, the large-scale correlations that imply 
a much larger particle horizon than the Big Bang sce- 
nario provides without inflation, and the nearly-scale- 
invariant Gaussian fluctuations. Although inflation was 
developed to explain these known properties of the Uni- 
verse, which are now probed wi th high precision by 
recent cosmol ogical observations (IKomatsu et alj l2011t 
Dunklev et al.ll2011b IKeisler et al.ll2011b lAnderson et all 
2012b iHicken et all 120091 : iKessler et al.ll2009t: IRozo et al.l 
2010ft . the model also has a new feature: the early, expo- 



nential expansion of space generates a stochastic back- 
ground of gravitational waves. In the near term, po- 
larization measurements of the cosmic microwave back- 
ground (CMB) present the most promising approach 
to detect these gravitational waves, which cause an 
odd-parity (B-mode) polari zation pattern on angular 
scales larger than a degre e (jSeliak fe Zaldarriagal 119971 : 
iKamionkowski et al.lll997|) . Detection (or non-detection) 
of these patterns will place strong constraints on the in- 
flation paradigm. 

In the slow-rol l appr oximation (for a review 
see ILiddle fc Lvthl 120001 ). the B-mode intensity 
is parametrized by the tensor-to-scalar ratio r, 
which is related to the energy scale V of infla- 
tion by V ~ (r/0.01) 1 / 4 x 10 16 GeV. For many 
classes of infla tionary models, r can be as large as 
0-01 < r < 0.1 (IBovle et all 12006ft . 

A combination of CMB-temperature-anisotropy 
measurements, baryon-acoustic-oscillation data, and 
supernova observations has given the most stringent 
limit to date, r < 0.2 at 95% confiden ce level (C.L.) 
nearl y limited by cosm i c variance (IKomatsu et al.1 
120111: IKeisler et all 120111: iDunklev et al.1 120111 ). In 
order to improve on these constraints significantly, 
direct observations of CMB polarization are re- 
quired. Thus far the best limit from C MB polariza- 
tion alone is r < 0.72 at 95% C.L. (jChiang et al.l 
12010ft . while many experiment s have ob served 
even-parity patterns (E- modes) dLeitch et al.l 
Montrov et alJl2006t ISievers et al.ll2007t IWu et"atl 



Bischoff et aTT 12008b iBrown et al l 120091 : lLarson et al.1 
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The Q/U Imaging ExperimenT (QUIET) observed 
the CMB from the ground between 2008 October 
and 2010 December. The observation site was the 
Chajnantor plateau at an altitude of 5080 m in the 
Atacama Desert in Chile. Two different receivers 
were employed, corresponding to center frequencies 
of 43 (Q-band) and 95 GHz (W-band). The results 
of the 43-GHz measurements hav e been published in 
IQUIET Collaboration et al.l (|2011f ) and included a mea- 
surement of the E'-mode power spectrum between I = 25 
and 475 and an upper limit on the tensor-to-scalar ratio 
of r < 2.1 at 95% C.L. In this paper, we report mea- 
surements of the CMB polarization power spectra for 
the 95-GHz data. We note that this experiment played 
the role of a pathfinder, demonstrating that monolithic- 
microwave-integrated-circuit (MMIC) arrays are capable 



of controlling systematic errors and achieving the sensi- 
tivity required to reach r < 0.01. 

QUIET was led by Bruce Winstein, who died in 2011 
February soon after observations were completed. His 
intellectual and scientific guidance was crucial to the ex- 
periment's success. 

2. INSTRUMENT 

In this section, we summarize the salient features of 
the 95-GHz instrum ent. For further details, we refer 
to separate papers (IQUIET Collaboration et al.l 12011b 
QUIET Collaboration 2012ft. hereafter referred to as 
QUIET 20111 and IQUIET 20121 respectively. Additional 
informat i on on the QU I ET in s trument is pr o vided in 
IBischoffl (12010ft: IBriziusI (I20H: IClearvl (l2010ft: IKusakal 
(12010ft: IMonsalvd (120101 ): INewburghl (12010L 12012ft : and 
iReevesI (|2012l ). 

The QUIET telescope consists of a 1.4-m side-fed clas- 
sical Drag onian antenn a that satisfies the Mizuguchi 
condition (|QUIET 2012ft . The Cosmic Background Im- 
ager (CBI) telescope mount was reused for the QUIET 
project. It provides three-axis motion: azimuth, eleva- 
tion, and rotation about th e optical axis, called "deck" 
rotation (|Padin et all 12002ft . The 95-GHz receiver com- 
prises 84 polarization-sensitive radiometers and six ra- 
diometers with differential-temperature sensitivity. The 
array sensitivity is 87/iK v / s to the CMB polarization. 
The instantaneous angular resolution is 11. '7 in full width 
at half maximum (FWHM) . The telescope field of view 
is roughly circular with a diameter of ~ 8° . 

The coherent QUIET radi ometers directly measure th e 
Stokes Q and U parameters (jQUIET 201 lb IQUIET 2012ft . 
The intensity, /, is also recorded by the same radiome- 
ters, but with significantly higher noise. One of the 
strengths of the QUIET design is excellent immunity to 
both 1/f noise from gain fluctuations and instrumen- 
tal spurious polarization (hereafter I-to-Q/U leakage). 
The median 1/f knee frequency of the radiometers is 
lOmHz, significantly below the typical scan frequency of 
45-100mHz, resulting in a negligible 1// noise contribu- 
tion. The fractional I-to-Q/U leakages are 0.2% for the 
monopole component, 0.4% for the dip ole component , 
and 0.2% for the quadrupole component (jQUIET 2012ft . 

The receiver and telescope mirrors are surrounded by 
an absorptive ground screen, eliminating major contri- 
butions from the 300-K ground emission. The upper 
component of the ground screen was installed in 2010 
January and eliminat ed two localiz ed far sidelobes with 
intensities 60 dB (jQUIET 20121) . which existed dur- 
ing the first few months of operation (from 2009 August 
through 2010 January). For the data from the early part 
of the season, we reject the part where the Sun entered 
either of these sidelobes. Scan-synchronous signal due 
to groun d emission is projected out of the maps in the 
analysis (IQUIET 2011ft . Possible remaining effects arc 
estimated as a systematic error (Section I5.3[) . 

3. OBSERVATIONS 

With the 95-GHz receiver, we observed from 2009 
August 12 until 2010 December 22 and accumulated 
7426 hours of dateP^I. Of these data, 72% were spent 

26 The instrument was in the nominal CMB observing con- 
figuration only between 2009 August 15 and 2010 December 17. 
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TABLE 1 

Data-selection Summary 





Observed Time 


Data Percentage 


Field 


(Hours) 


ML 


PCL Both 


CMB-1 


1855 


69.7 


64.4 57.7 


CMB-2 


1444 


73.1 


67.1 61.2 


CMB-3 


1389 


64.4 


58.8 52.6 


CMB-4 


650 


72.1 


65.4 60.4 


Total 


5 337 


69.5 


63.5 57.6 


Note. 


Fraction 


of data 


selected for 



each field by each pipeline. The last column 
shows the fraction simultaneously selected by 
both pipelines. 

on CMB observations, 14% on Galactic fieldf£], 13% on 
calibration sources, and 1% on incomplete observations 
due to obvious instrumental problems such as a lack of 
telescope motion. We observed 24 hours per day, except 
for interruptions due to a variety of factors such as high 
wind, heavy snow, power outages, and instrumental fail- 
ures. Our full-season operating efficiency was 63%. For 
the CMB measurements, we selected fo ur low-foregro und 
sky fields, denoted CMB-1, 2, 3 and 4 (IQUIET 20111 ). In 
total, we collected 5337 hours of CMB data with the 95- 
GHz receiver (Tabled]). 

Each observation consists of a series of constant- 
elevation scans, hereafter collectively called a CES. The 
scans are in the azimuth direction with a half amplitude 
of 7.5° on the sky. Diurnal motion of the sky causes the 
field to drift through the field of view. After the target 
has drifted 15° on the sky, we adjust the azimuth and 
elevation to retrack the field and begin a new CES. Each 
individual CES thus scans over an area of ~ 15° x 15°. 
Due to the field of view of ~ 8° and the fact that the sky 
does not always drift orthogonal to the scan direction, 
a larger area is observed in practice. The deck angle 
is changed by 45° each week, providing a large degree 
of immunity to spurious -B-modes induced by I-to-Q/U 
leakage. 

4. CALIBRATION 

The instrument calibration procedure for the 95-GHz 
obser vations is similar to that used for the 43-GHz 
data (jQUIET 20111: IQUIET 20121) . The instantaneous 
beam point-spread function is derived from observations 
of Taurus A (hereafter Tau A). The resulting beam func- 
tion has a width o f 11. '7 FWHM with a small non- 
Gaussian correction (| QUIET 201 2D . The telescope point- 
ing model is calibrated with a set of astronomical objects: 
Tau A, Jupiter, RCW 38, the Moon, and the Galactic 
center. The residual random scatter after applying all 
pointing corrections is FWHM. To correct for this, 
we convolve the beam window function with the residual- 
pointing-scatter term, and obtain an effective point- 
spread function of 12f8 FWHM. The detector angles (i.e., 
the orientations of the polarization responses) are cali- 
brated to 0?5 precision with the combination of Tau A ob- 
servations for absolute- angle determination and a sparse- 
Different configurations were used between 2009 August 12 and 15 
and between 2010 December 17 and 22 to calibrate and character- 
ize the instrument. 

27 The ana lysis of the Galactic observations is in progress (see 
I Wehusl l|2012l ) for preliminary maps), and final results will appear 
in a future publication. 



wire-grid calibrator (jTaiima et all [20121: IQUIET 20121) 
for relative angle determination. The considerable im- 
provement in the detector angl e precision co mpared to 
the previous 43-GHz analysis (IQUIET 20111) is due to 
a mo re accurate catalog value of Tau A (jAumont et al.1 
E0U as well as an improved wire grid calibration. 
Large and small sky dips (elevation nods of ±20° and 
±3° amplitudes, respectively) modulate loading from at- 
mospheric emission and allow us to measure the frac- 
tional I-to-Q/U monopole leakage with 0.3% precision 
per calibration, while Jupiter measurements are used 
to measure the higher-order leakage terms (i.e., dipolc 
and quadrupole) and to confirm the sky-dip monopolc- 
lcakagc results. The detector responsivities are cali- 
brated using Tau A and sky-dip data as well as the mea- 
surement using the sparse wire grid. The typical rcspon- 
sivity is found to be 3.1 mVK -1 in antenna temperature 
units. 

5. DATA ANALYSIS 

The analysis procedure used for the 95-GHz 
dat a reduction follows closely the 43-GHz analy- 
sis (jQUIET 20111 ) , and we refer the reader to this publica- 
tion as well as recent Ph. D. theses (iBuder 2012: C hinonel 
l2lHltlDumouhnll2l^ for full 

details. We have implemented two independent analy- 
sis pipelines, one based on a maximum-likelihood (ML) 
technique and the other on a pseudo-Cf (PCL) cross- 
correlation technique. The most important improve- 
ments since the previous publication are, for the ML 
pipeline, an adaptive filter procedure in which the fil- 
ter parameters depend on the data quality of the specific 
data segment, as well as a pseudo-C^ null-test estima- 
tor, allowing for many more null tests; and, for the PCL 
pipeline, a different and more robust data division for 
the cross-correlatiorEl, taking further advantage of the 
scanning strategy. 

The process to extract cosmological results from raw 
time-ordered data (TOD), containing measurements of 
the Stokes Q and U parameters as well as the tele- 
scope pointing information, can be summarized in three 
steps: TOD pre-processing, map making, and power- 
spectrum and parameter estimation. The TOD pre- 
processing involves estimating and applying calibration 
factors, characterizing the detector noise, and applying 
high-pass, low-pass and azimuth filters to minimize the 
effects of atmospheric fluctuations, far sidelobes, excess 
high-frequency instrumental noise, and ground pickup. 
Then, sky maps are generated by projecting the Q and U 
intensities into Galactic coordinates, taking into account 
the telescope pointi ng informatio n, using standard map- 
making equations (jQUIET 20111 ). Figure Q] shows the 
maximum-likelihood Stokes Q and U maps of the CMB- 
1 field generated with the ML pipeline. Power-spectrum 
estimation is performed with one of two techniques, de- 
pending on the pipeline. The ML pipeline implements 
a standard Newton -Raphson maximum-likelihood solver 
(jBond et al.lll998l ). while the PCL p ipeline implements 
the MASTER pseudo- fi algorithm (jHivon et all l200l 
lHansen &; Gor ski 2003). Prior to power-spectrum esti- 

28 We cross-correlate among 40 subsets of data. Each subset 
corresponds to a specific boresight azimuth and deck range. There 
are five azimuth and eight deck ranges. 
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mation, both pipelines mask Ccntaurus A, and the PCL 
pipeline also masks Pictor A. 

In the following we describe the data selection, analysis 
validation and systematic-error assessment. In optimiz- 
ing the analysis configuration, it is important that the op- 
timization process itself does not introduce experimenter 
biases, for instance by removing purely statistical fluc- 
tuations in the data selection. QUIET is the first CMB 
experiment to have adopt ed a strict blind-analysis policy 
(Kl ein fc Roodmanir2005[ ). in which all data-selection cri- 
teria, filters, and calibrations are adjusted and finalized, 
and the systematic errors are assessed prior to looking 
at any cosmological power spectru m. This process was 
described in detail in lQUIET 20111 and we have adopted 
the same policy for the 95-GHz analysis. 

5.1. Data Selection 

Each QUIET radiometer provides four output channels 
("detector diodes"), resulting in a total of 336 output 
channels from 84 polarization-sens itive radiomet ers. We 
use 308 good channels for analysis ([QUIET 2012[ ). Start- 
ing from the resulting data for all CESes, we define two 
different classes of data-selection criteria. In the first 
class, we impose criteria that select or reject an entire 
CES. These include the criteria based on atmospheric 
conditions, instrument malfunctions, or unusual condi- 
tions for the temperature regulation in the focal plane. 
In the second class, we apply selection criteria to individ- 
ual detector diodes in each CES (CES-diodcs). For in- 
stance, a CES-diodc is rejected if: 1) the measured noise 
properties show poor agreement with the noise model; 
2) the 1// knee frequency is anomalously high; 3) the 
white-noise level is non-stationary; 4) there are glitches 
in the time domain or strong spikes in the Fourier do- 
main; or 5) there is evidence of a large scan-synchronous 
signal. 

Table Q] lists the fractions of data that satisfy the cri- 
teria and are used for map making and power-spectrum 
estimation. 

5.2. Analysis Validation 

Having defined our data-selection criteria and filters, 
we need to validate the accepted data set and analysis pa- 
rametei'fQ O ur most valuab le tool for this is a so-called 
null-test suite (jQUIET 20 111 ). In each null test, the full 
data are split into two subsets. From these, we make indi- 
vidual sky maps, mi and m2, as well as the corresponding 
difference map, mdifr = {m\ — n%2)/2. By design, the true 
sky signal cancels in this map, and the result should be 
consistent with noise. We therefore compute the EE and 
BB power spectra of this map, and check for consistency 
with the zero-signal hypothesis by comparing to simula- 
tions. In the current analysis, the null suite consists of 
32 and 23 tests for the PCL and ML pipelines, respec- 
tively, with each test targeting a possible source of signal 
contamination or miscalibration. These are selected to 
be highly independent; a statistical correlation between 
null power spectra of two different null-test divisions is 
typically 0.05. 

29 Note that the data-selection criteria are improved through 
an iterative process of applying the analysis- validation metrics. In 
the current analysis, ~ 50 different configurations were considered 
before reaching the final configuration. 



TABLE 2 
Validation-test Summary 



Field 


Mean of Xnull 
ML PCL 


Total 
ML 


^null 

PCL 


^null 

ML 


Outlier 
PCL 


CMB-1 


0.40 


0.26 


0.31 


0.78 


0.63 


0.14 


CMB-2 


0.54 


0.46 


0.08 


0.06 


0.02 


0.40 


CMB-3 


0.42 


0.74 


0.31 


0.50 


0.72 


0.45 


CMB-4 


0.06 


0.08 


0.43 


0.76 


0.80 


0.18 



Note. — Results of the three predefined validation 
tests using the mean of Xnull i the sum of X n ip an d the 
worst outlier of X nu ii- All values are PTEs defined such 
that a large deviation from zero results in a small PTE. 



For each power-spectrum bin b, we calculate the statis- 
tic XmulO) = C£ lull /cr fc , where C fc nu11 is the observed dif- 
ference power spectrum and en, is a Monte Carlo (MC) 
based estimate of the corresponding standard deviation. 
We evaluate both Xnuii and its square for all b; x nu ii is 
sensitive to systematic biases in the null spectra, while 
X nu n is more responsive to outliers. 

Prior to the analysis, we defined three critical tests that 
had to be passed before continuing to cosmological anal- 
ysis, based on 1) the mean value of Xnuii, 2) the sum of 
Xnuii' ana - 3) the maximum of Xnuii' au computed includ- 
ing the entire suite of EE and BB null power spectra. 
A given analysis configuration passes when these statis- 
tics are consistent with the null hypothesis. Table [5] lists 
the probabilities to exceed (PTE) for the final configura- 
tion, and Figure [2] shows the PTE distribution of X nu ii- 
The PTEs are defined such that a large deviation from 
zero results in a low PTE. This corresponds to two-sided 
PTEs for the mean of Xnuii and one-sided PTEs for the 
total Xnuii an( i the- Xnuii outlier. The mean of the Xnuii 
distributions over all fields is -0.018 ±0.015 for the PCL 
pipeline and 0.003 ± 0.017 for the ML pipeline. We do 
not detect any bias with our final analysis configuration. 

We also generate 1000 random null divisions and com- 
pare the widths of the resulting Xnuii distributions be- 
tween data and MCs using the PCL pipeline. We find 
these to be consistent, and we verify our estimate of the 
statistical uncertainty in each multipole bin with a pre- 
cision of 3.1%. Finally, we evaluate the differences of 
non-null spectra among the fields, before looking at indi- 
vidual non-null spectra. These differences are consistent 
with the hypothesis of statistical isotropy (i.e., each field 
has the same underlying power spectrum), with a PTE 
of 0.15. 

5.3. Systematic Errors 

We study the contributions from in strumental sy stem- 
atic errors using the methodology of IQUIET 201ll The 
main effects considered are 1) uncertainties in absolute 
responsivity and the window function, 2) I-to-Q/U leak- 
age, 3) uncertainties in polarization angles, relative re- 
sponsivities, and pointing, and 4) residual contamination 
from scan-synchronous signals and far sidelobes. In each 
case, we set up an empirical model of the systematic ef- 
fect and propagate this through the PCL pipeline. The 
results from these calculations are summarized in Figure 

El 

The most important conclusion is that the systematic 
errors in the BB spectrum are very small. For the multi- 
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Fig. 1. — QUIET CMB polarization maps of the CMB-1 field in Galactic coordinates at 95 GHz. The left (right) panel shows Stokes Q 
(U), where the polarization angle is defined with respect to the Galactic North Pole. Note the coherent vertical/horizontal patterns in the 
Q map, and the diagonal patterns in the U map; these are the expected signature of a pure i?-mode signal. No filtering has been applied 
to this map beyond subtracting the very largest angular scales (I > 25), to which QUIET is not sensitive. 

2 , 

pole range relevant for estimation of the tensor-to-scalar 
ratio, £ ~ 100, each effect is smaller than or compara- 
ble to the signal corresponding to r ~ 0.01, the lowest 
level ever reported in the literature. It is also noteworthy 
that this limit improves on that reported for the 43-GHz 
data (r < 0.1; QUIET 2011) by an order of magnitude. 
This is due to improved rejection of I-to-Q/U leakage, 
better detector- angle calibration, and lower levels of side- 
lobe contamination resulting from the installation of the 
upper parts of the ground screen. 

For the EE power spectrum, the systematic error bud- 
get is dominated by uncertainties in the multiplicative 
responsivity calibration. The total uncertainty is 8%, al- 
most equally contributed from three dominant sources: 
the uncertainty of t he polarization flux of Tau A (5%; 
iWeiland et al J 120 111) , the uncertainty in the beam solid 
angle (5%), and the uncertainty associated with model- 
ing the time variation and relative responsivity among 
the detector channels (4%). This translates into an un- 
certainty of 17% in the power spectrum. For compar- 
ison, the statistical uncertainty in the EE spectrum is 
about 8% of the central value at its minimum around 
I ~ 400. It is important to note that the responsivity ef- 
fect is purely multiplicative and therefore cannot create 
spurious B-mode signal. The uncertainty of the window 
function is another multiplicative factor highly correlated 
among different £ bins, with the magnitude dependent on 
I. The uncertainty comes from both the beam window 
function and the smearing factor due to the pointing er- 
ror, and is listed in Table [3l These errors are smaller 
than the EE statistical uncertainties. 

The dominant systematic uncertainty for EB is due 
to calibration errors in the detector polarization an- 
gles. To first approximation, an error in the absolutc- 
polarization-angle calibration of dip induces a spurious 
EB spectrum proportional to ~ Cf E sin 28i/j, and a BB 
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Null-test PTE 

Fig. 2.— Null-test-PTE distributions for y 2 „ for both the ML 

^-null 

and PCL pipelines. Each is consistent with the uniform expecta- 
tion. 

a similar manner. The calculations summarized in Fig- 
ure [3] capture both these effects through simulations. As 
seen in this figure, these polarization-angle uncertainties 
lead to systematic errors almost as large as the statis- 
tical errors for EB around £ ~ 400, while for BB they 
are small everywhere and comparable to other sources of 
systematic errors. Table[3]lists the total systematic error 
for the EB power spectrum. 



6. 



POWER SPECTRA AND COSMOLOGICAL 
PARAMETERS 



spectrum proportional to ~ Cf sin 28ip. Uncertainties 
in the relative polarization angles among detectors con- 
tribute to the systematic errors in EB and BB spectra in 



The measurements of the EE, EB and BB power spec- 
tra are tabulated in Tablc[3j and plotted in FigurclU The 
EE spectrum is strongly signal-dominated up to £ ~ 800, 
and three acoustic peaks are clearly traced. Both the 
BB and EB spectra are consistent with zero within the 
estimated statistical and systematic uncertainties. The 
dominant EE power is also visible in the maps shown in 
Figure [T] Note that these maps have not been filtered, 
except by subtracting the very largest scales (£ < 25), 
to which QUIET is not sensitive. One can see a dis- 
tinct vertical-horizontal coherent pattern on small angu- 
lar scales in the Stokes Q map, and a similar diagonal 
pattern in the U map. This is the expected signature of 
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Fig. 3. — Summary of systematic error assessment for EE (left), BB (middle), and EB (right). The red bars indicate the statistical 
uncertainties in each bin. Blue, green, and purple points correspond to three categories of systematic errors: I-to-Q/U leakage; polarization 
angles (absolute and relative), relative responsivities and pointing error; and the residual scan-synchronous signals and far sidelobes. The 
gray band along the ACDM curve in EE corresponds to the uncertainties of multiplicative factors: absolute responsivity and the window 
function. For BB, all systematic errors are below the level of r ~ 0.01 at i ~ 100. For EE the dominant systematic error is uncertainty 
in the absolute responsivity, which is a purely multiplicative effect. For EB, the dominant systematic is caused by uncertainties in the 
polarization detector angle. 

an E-modc signal. 

The results from the two pipelines are consistent with 
each other. The most noticeable difference is a single 
overall multiplicative factor, which is only relevant in 
evaluating the consistency of the EE power spectra. This 
factor comes from different responsivity modeling and is 
consistent with the systematic error budget discussed in 
Section 15.31 

When assessing the consistency of the EE power spec- 
trum with the ACDM prediction, it is convenient to fac- 
tor the spectrum measurement into an overall amplitude 
and the spectral shape of the acoustic peaks. We fit 
a free amplitude, q, relative to the EE spectrum pre- 
dicted by the best-fit seven- year WMAP ACDM param- 
eters (jKomatsu et al.ll2011[ ) to the spectrum from each 



pipeline, and find q = 1.22 ± 0.04(stat)t ;i7(syst) and 
q = 1.35 ± 0.05(stat)±gJ|(syst) for the PCL and ML 
pipelines, respectively. These values are consistent with 
the ACDM prediction of q = 1, and correspond to PTEs 
of 0.20 and 0.06, respectively. Figure [5] provides a spec- 
tral shape comparison. Here we see that the measured 
EE spectrum rescaled to q = 1 accurately traces the 
first three acoustic peaks predicted by the ACDM model. 

We assess the overall consistency with the ACDM 
hypothesis by calculating a total \ 2 relative to ACDM 
(and relative to Cf B = Cf B = 0), taking into account 
the systematic uncertainties due to the responsivity cali- 
bration in EE and the systematic error in EB primarily 
due to detector polarization angles. The former is in- 
corporated by introducing a nuisance parameter for the 
absolute responsivity constrained by a Gaussian distri- 
bution with a standard deviation equal to the assigned 
systematic error. The latter is incorporated by mod- 
eling the EB systematic error as sC EB ' syst , where the 
scale factor s is constrained by a Gaussian with a = 1 
and (j^ B ^y &t i s the systematic error estimated in Sec- 



tion [573J this means we assume the systematic errors are 
completely correlated among different £ bins. The sys- 
tematic errors in BB are negligibly small. Including the 
systematic error contributions, we find x 2 °f 67.3 and 
70.2 for the PCL and ML pipelines, respectively With 
57 degrees of freedom, these values correspond to PTEs 
of 0.16 and 0.11, respectively; the derived spectra are 
consistent with ACDM. 

Since we find no significant excess in the BB power 
spectrum, we place an upper limit on possible BB power 
in each bin. The ML pipeline calculates the upper limit 
by the 95% integral of the positive part of the likeli- 
hood, while the PCL pipeline adopts a frequentist-based 
hypothesis-testing method. Specifically, the upper limit 

' 10), 



obs 



fi is defined by 0.05 

where p(- ■ ■ \fi) and p{- ■ ■ |0) represent p- values of the null 
hypothesis with power /x and an alternative hypothesis 
with zero power, respectively. The parameters fi, q^ and 
g° bs correspond to the bandpower Cb, the test statistic 
for upper limit defined in iCowan et al.l ()2011[ ) , and the 
test statistic q^ calculated for the observed bandpower 

Cb, respectively. Table |3] lists the derived upper limits. 

We constrain the tensor-to-scalar ratio r using stan- 
dard likelihood methods and including only the BB spec- 
trum at low multipoles (26 < I < 175). For sim- 
plicity, we consider only the amplitude of a BB tem- 
plate computed with the standard ACDM concordance 
pa rameters, and fix the tensor spe ctral index to n t = 
(jQUIET 20111 : iChiang etaD 120101 ) . In constraining r, 
the uncertainty of the responsivity calibration is elimi- 
nated by simultaneously fitting EE and BB power spec- 
tra using the ACDM templates. We define the fit func- 
tion as Cf E {q) = qCf Efid and Cf B {q,r) =rqCf Bm . 

Here C EEfid and cf Bfid denote the fiducial ACDM SB- 
power-spectrum template and the Bi?-power-spectrum 
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Fig. 4. — The QUIET 95-GHz power spectra, co-added over all 
four CMB fields. The panels show the EE (top), BB (middle) 
and EB (bottom) spectra, and the insets show the low-£ region in 
detail. The central I values from the two pipelines arc slightly offset 
for display purposes. N ote t hat the error bars indicate statistical 
errors only; see Section I5.3I for a discussion of systematic errors. 
Typical correlations among neighboring bins are ~ —0.1. The full 
set of three spectra are consistent with the EE spectrum predicted 
by the ACDM model and Cf B = Cf B = 0. 
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Fig. 5. — Comparison of the QUIET EE spectrum with current 
best-fit ACDM model after scaling the absolute responsivity to 
q = 1 . The central I values from the two pipelines are slightly offset 
for display purposes. Typical correlations among neighboring bins 
are ~ —0.1. The results from the two pipelines are consistent, and 
the shape of the QUIET EE spectrum is in excellent agreement 
with the ACDM model. 

template with r = 1 and n t = 0, respectively. Note that 
this does not imply that we use EE power to constrain 
the tensor modes, as the ACDM EE template only con- 
tains the scalar contribution. This method exploits the 
fact that r is by definition a ratio and does not depend 
on the common overall scaling factor. From the simulta- 



neous fit, the ML pipeline finds r = 1.1 



+ 0.0 



with a 95% 



C.L. upper limit of r < 2.8, and the PCL pipeline finds 
r = 1.2+Qg, corresponding to an upper limit of r < 2.7. 
The systematic uncertainty is negligible, at the level of 
r = 0.01. 

7. FOREGROUNDS 

We assess the level of diffuse foregrounds, in particular 
synchrotron radiation and dust emission, as additional 
sources of systematic errors. Contamination from resid- 
ual point sources is negligible. An estimate using the 
point-source component of the Planck Sky Model (PSM; 
iDelabrouille et all 120121 PSM vl.7.4) yields a limit of 
Cg < 1.4 x 10~ 6 /jK 2 over the entire I range without 
masking any sourc es^]. An estimate based on a source- 
population model ([Tucci k, ToffolattH |2012[ ) relative to 
our nominal point-source mask results in an even lower 
level, Cg ~ 5 x 10~ 7 ^K 2 . Both are well below our sta- 
tistical uncertainty. 

Considering synchrotron radiation, we note that the 
43-GHz QUIET observations have already resulted in 
strong constraints on any syn chrotron comp onent in each 
of the QUIET CMB fields (| QUIET 20111 ). Except for 
the single case of the EE spectrum at I < 75 mea- 
sured in CMB-1, no evidence of any contamination was 
found. These results allow us to constrain any con- 
tribution from synchrotron emission at 95 GHz by ex- 
trapolation. Adopti ng a spectral index of /3 S = —2.7 
(|Dunklev et al.ll2009l) . we estimate the EE (BB) excess 
power to be 0.011 ± 0.003 fiK 2 (0.001 ± 0.002 ^K 2 ) for 
the first bin of the CMB-1 spectrum, which is negligible 
compared to statistical errors. 

In order to constrain contamination from dust emis- 

30 Note that this limit is given in units of Cg, not Cg l)/2n. 
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sion, wc adopt the thermal-dust component of the PSM 
as a template; the PSM predicts that other sources of 
contamination are subdominant at 95 GHz in the QUIET 
fields. We estimate the dust power contribution in our 
fields by evaluating both the PSM power spectrum and 
the PSM-QUIET cross-spectrum using the PCL pipeline. 
The possible contamination is only relevant in the first 
bin (25 < I < 75) of the field CMB-1. In this bin, the 
PSM power amplitude is 0.087 ^K 2 (0.070 /zK 2 ) for the 
EE (BB) spectrum, while the corresponding cross power 
is 0.060 ± 0.035 /zK 2 (0.016 ± 0.027/zK 2 ). Taking into 
account the relative weights of the individual fields, we 
therefore estimate that the dust-emission contribution to 
the first EE bin in the final co-added spectrum (Table 
[3]) is < 0.04 /iK 2 , more than a factor two smaller than 
the statistical uncertainty. All other spectra and mul- 
tipolc ranges have negligible contributions. Fitting the 
PSM model as a template to CMB-1 in the map do- 
main using the ML pipeline, we find a best-fit amplitude 
of A — 0.62 ± 0.21. This corresponds to a 3cr corre- 
lation with the thermal-dust PSM component, which at 
the same time agrees with the PSM prediction (A = 1) at 
1.8(7. Consistent results are obtained by taking the ratio 
of the cross-power to the PSM power including the full 
multipolc range, with an amplitude of A = 0.66 ± 0.25. 
The three other fields all have best-fit amplitudes consis- 
tent with zero. We note as a caveat that the uncertainty 
in the PSM itself is not taken into account in this analy- 
sis, and the results depend critically on this model as the 
detected foreground levels are well below the statistical 
errors of the measured power spectra themselves. 

8. CONCLUSIONS 

We have presented the CMB polarization power spec- 
tra from the 95-GHz QUIET observations. The EE spec- 
trum has been measured between i = 25 and 975, and 
the first three acoustic peaks were seen with high signal- 
to-noise ratio, consistent with ACDM predictions. The 
BB spectrum was found to be consistent with zero, with 
a 95% C.L. upper limit on the tensor-to-scalar ratio of 
r < 2.7 (PCL) or 2.8 (ML), depending on pipeline. In 
Figure [6l we provide an up-to-date overview of the cur- 
rent state of the CMB polarization field, comparing the 
results from various experiment J^l. In one of the fields, 
we found a correlation with the dust component of the 
Planck Sky Model. The excess power due to this com- 
ponent was still small compared to the statistical errors 
of the power spectra. Finally, we have demonstrated the 
lowest level of instrumental systematic errors to date. We 
conclude by noting that part of the role of this experi- 
ment was to serve as a pathfinder to demonstrate that 
MMIC arrays were capable of reaching r < 0.01; this has 
been successfully achieved. 
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QUIET Polarization Power Spectra 
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±0.03 


176 


225 


0.53l°;g (1.09) 


„ ^ ~-4-0 99 
u - lo -0.22 


0.39i° j2 (0.95) 


17+ - 24 

u ' 1 '-0.23 


±0.02 


226 


275 


-0.59±°;|J (0.52) 


~ . , 4_fl "i^ 
nil -fu.dj 

U - li -0.34 


-O.6OI031 (0.49) 


12+ 033 
u ' lz -0.32 


±0.08 


276 


325 


-0.4+°;« (1.0) 


^ -. -1-0 fi 


-0.2+_1i (1.0) 


1+ 06 


±0.29 


326 


375 


-0.5l°;? (1.2) 


_l_0. 8 


-o.ilg:? (i.5) 


2.o+g:i 


±0.59 


376 


425 


-0.lll;° (2.1) 




0-411.0 (2.5) 


0.911.1 


±0.75 


426 


475 


1.6lJ;S (4.4) 




1.211J (3-7) 


2.5111 


±0.55 


<± i 






U - 9 -1.6 


-2.5ll'7 (2.3) 


2 +1 + 
u - z -1.3 


-1-0 '30 


526 


575 


4.21 2 ;?, (9 . 2) 


3+ 20 


2.3ll; 4 3 (6-6) 


-i.6i^ 9 o 


±0.28 


576 


625 


3-8±|;5 (10.5) 


o 4 +2.7 


-l.ll 3 2 ;i (5.6) 


i.8ll:l 


±0.70 


626 


675 


(18.3) 




6.31+1 (13-8) 


4 0+3.5 


±1.25 


676 


725 


3t£ (16) 


ill 


8l 6 6 (18) 


6 7+ 4 - 8 


±1.40 


726 


775 


-8±1 (13) 


9 + 6 


ll 8 s (17) 


_11+5 
lL ~6 


±1.06 


776 


825 


7li| (31) 


-19±g 


511J (26) 


-1411 


±0.68 


826 


875 


10i^ (44) 


15+12 


2ll 5 4 (31) 


_q+l0 


±0.55 


876 


925 


46+^ (93) 


4 +19 

4 -19 


-3111 (37) 


Q+15 


±0.93 


926 


975 


-52lg (44) 


n +25 
u -25 


-28±aa ( 41 ) 


,-+19 
11 -19 


±1.59 



Note. — Tabulated values are given in CMB thermodynamic units of /iK 2 , scaled as C{£(£ + l)/27r. 
We present the results from both the ML and PCL pipelines; they are in excellent statistical agreement. 
The column EE/q shows the EE power spectrum normalized to q = 1, as plotted in Figure [S] The fit 
value of q is also shown in the table, where the first and second errors are statistical and systematic, 
respectively. The column of BB-power-spectrum values also provides 95% confidence level upper limits 
in parentheses. We also list two relevant systematic-error contributions besides the uncertainty from 
the responsivity calibration: the fractional error due to the uncertainty of the beam window function, 
and the total systematic error in the EB power spectrum in units of fiK 2 . Note that they are both 
highly correlated among t bins. 



